Stochastic Processes with Balking in the 
Theory of Telephone Traffic* 

By LAJOS TAKACS 

(Manuscript received January 16, 1961) 

// is supposed that at a telephone exchange calls are arriving according 
to a recurrent process. If an incoming call finds exactly j lines busy then 
it either realizes a connection with probability p } or balks with probability 
Qi ("Pi + Qj " !)• The holding times are mutually independent random 
variables with common exponential distribution. In this paper the stochastic 
behavior of the fluctuation of the number of the busy lines is studied. 

I. INTRODUCTION 

Many results in telephone traffic theory (and elsewhere) may be 
unified by the introduction of balking. A call is said to balk if for any 
reason it refuses service on arrival. A mathematical model for balking is 
constructed by assigning a probability to balking dependent only on the 
state of the system; if an incoming call finds exactly j lines busy, then 
it realizes a connection with probability pj and balks with probability 
Qj (Pi + (jj = 1). Thus if Pj = 1 (j = 0, 1, ■ ■ •) the system is one with 
an infinite number of lines and with no loss and no delay, the ideal for 
for any service, while if p } ■ = 1 (j = 0, 1, ■ • •, m - 1) and p> = 
(j = m , m + 1, ■ • • ) the system is a loss system with m lines. 

This balking model is examined here for recurrent input and expo- 
nential distribution of holding times. More specifically, the call arrival 
times are taken as the instants n , r-> , • • •, r„ , • • •, where the inter- 
arrival times 6 n = r«+i - r„ (n — 0, 1, • ■ •; t = 0) are identically 
distributed, mutually independent, positive random variables with dis- 
tribution function 

(1) 



* Dedicated to the memory of my professor Charles Jordan (December 16, 
1871-December 24, 1959; 
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The holding times are identically distributed, mutually independent 
random variables with distribution function 

H(x) -\0 if*<0. (2) 

The holding times are independent of the {r„( as well. 

Let us denote by £(t) the number of busy lines at the instant t. De- 
fine h n = £ (t„ — 0); that is, £„ is the number of busy lines immediately 
before the arrival of the nth call. The system is said to be in state E k 
at the instant t if £(/) = k. Let us denote by m the smallest nonnegative 
integer such that p,„ = 0. If p } > (j = 0, 1, 2, • • ■ ) then m = oo. 

In the present paper we shall give a method to determine the dis- 
tribution of £„ for every n, the distribution of £(/) for finite / values, 
and the limiting distributions of £„ and £(t) as n — > °° and / — » <» re- 
spectively. Further, we shall determine the stochastic law of the transi- 
tions E k -+ E k+1 (k = 0, 1,2, •••)• 

II. NOTATION 

The Laplace-Stieltjes transform of the distribution function of the 
interarrival times will be denoted by 



vis) =E{<r s9 «} = r G -° x dF(x), 

Jo 



which is convergent if (R(s) ^ 0. The expectation of the interarrival 
times will be denoted by 

a = E{0„} = f xdF(x). 
Jo 

Let PU„ = ft) = P k (n) and P{£(0 = A-} = P k (t). Define 

U k (s) = [ e~ st P k (t)dt, 

Jo 

which is convergent if (R(s) > 0. Let 

lim P k 0,) = P k and lim P k (t) = P k *, 

n-»« J-»oo 

provided that the limits exist. 
Define 

C r -Il( r ^ n ) (r- 0,1,2,...), (3) 
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where the empty product means 1; that is, C = 1. We shall also use 
the abbreviation 

<p r = p(r/0 - f C" dF(x) (r = 0, 1, 2, • • • ) . (4) 

Jo 

Denote by M k (t) the expected number of calls occurring in the time 
interval (0,/) which find exactly h lines busy. The expected number of 
transitions E k — > E k+i occurring in the time interval (0,/] is clearly 
p k M k {t). Denote by N k (t) the expected number of transitions E k+l —> E k 
occurring in the time interval (0,t\. 

Let G k (x)(k = 0, 1, 2, • • • ) be the distribution function of the time 
differences between successive transitions Et-i — ► E k and E k — > E k+l , 
while R k (x)(k = 0, 1, 2, • • •) is the distribution function of the time 
differences between consecutive transitions Ek —* Ek+i . If £(0) =0 
then we say that a transition £"_i — > E takes place at time t = — 0. 
Write 



7*0) = f e-" dGk(x) 



'0 

and 

Pkis) = f° e-*dR k (x) 

Jo 

which are convergent if 6\(s) ^ 0. 

III. PREVIOUS RESULTS 

3.1 A. K. Erlang 

Erlang 1 has proved that, if { r„} forms a Poisson process of intensity 
X — that is, F(x) = 1 — c~ Xx for x ^ — and further, p> = 1 when j < m, 
Pj = when j ^ m, then 



(x/m: 



y- (X/MJ 

In this case P* = 1\* (ft = 0, 1, • • •, m) also holds. This is the simplest 
loss system. 
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3.2 Conny Palm 

Palm 2 has generalized the above result of Erlang for the case when 
{t,,} forms a recurrent process and otherwise every assumption remains 
unchanged. Palm has proved that 

P m = 5 /!a ; , (6) 



7m\l> 
t=i\r)Cr 



z 



where C\ is defined by (3). In this case the complete limiting distribu- 
tions [Pk] and {Pk*} have been determined by Pollaczek, 3 Cohen, 4 and 
the author. 6,6 The transient behavior of the sequence {£,} was determined 
by Pollaczek 3 and Benes, 7 and the transient behavior of the process 
{{(01 b y BeneS 8 and by the author. 9 

3.3 The Infinite Line Case 

The case when p } = 1 (j = 0, 1, 2, • • • ) has been investigated by the 
author, 10 '" who has proved that 

Pk = E (-l)'"* (fc) Cr (* - 0, 1, 2, • • • ) (7) 

and, if F(x) is not a lattice distribution and if a < <*>, then the limiting 
distribution {Pk*} exists and 

Pk* = P^ (h = 1,2, ■••), 
Kau 

(8) 



V = l--I^. 



The transient behavior of the process {£(/)} is also treated in Refs. 10 
and 11. 

3.4 The Case p = 1, Vi = p (j = 1, 2, • • 0, ft - <7 = 1, 2, • • •) 

This case, where p + q = 1, plays an important role in the theory of 
particle counters and has been investigated by the author, 1 "' who has 
found that 

pZ(-pYCr 

Po = =^ (9) 

1 -«E(-Pm 
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and 

E(-ir<(;:Vv 

p fc = r=* = — W— (fc = l,2, ...). do) 

r=() 

If F(z) is not a lattice distribution and if a < <x>, then the limiting 
distribution {P**j exists and 

P *H*-7TXTr- (fc = 0,1,2,.-.), 



/v = i - - E 



p*p* 



a^ A=n (A - + 1 ) 

The transient behavior of the process {£(/)} is a,so treated in Kef. 12. 

3.5 The Distribution Function G*(x) 

This function plays an important role in the investigation of overflow 
traffic. In the infinite line case, i.e., when pj = 1 (j = 0, 1, 2, • • •), 
Palm 2 has proved that y k (s)(k = 0, 1, 2, •••) satisfies the following 
recurrence formula : 

7*(«) = i ,„ N , /„ ■ n k* - If V"). (1 -) 

1 — 7*_i(s) + 7a-i(S + m) 

where y (s) = p(s). Palm has obtained yds) explicitly when \t„] is a 
Poisson process; that is, <p(s) = X. ( X + s). Then 

sis +m)-"[«+ <•/ - Om] 



s© 



/ \ _ >=o \.// V (1*3) 

The general solution of the recurrence formula (12) is 

*<•> = * + r;/iVw- L , r r fa + ^i An ^ = o, i, 2, . . . ) , (u) 

h\ r /fiL *(« + t/0 J 

where the empty product means 1. The formula (14) is proved in Refs. 
10 and 11. 
In the particular case p = 1, pj = p (j = 1, 2, •••)></> = 



(16) 
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q(j = 1,2, ■ • •), where p + q = 1, the Laplacc-Stieltjes transform y k (s) 
has been given explicitly in Ref. 12. We have 

Y*(s)=-^rv (k = 0,1,2,...), (15) 

where D (s) = 1 and 

dm = {p t (*) n p'fl+tl 

_ g[l - y?(s)] y« //A V ( _ 1 )> TJ f l - <PJS + Zm) "[ 

if fc = 1,2, •••. 

IV. THE TRANSIENT BEHAVIOR OF {£„} 

It is easy to see that the sequence of random variables (£„} forms a 
homogeneous Markov chain with transition probabilities 

p jh - P{$ n+ i = k ||„ = j) = f T jk (x) dF(x) , (17) 

•'0 

where 

Tjk(x) = 

Pi + l \ e~ k >"(l - e -n i+1 - k + qj ({) r**(l - e^)*"* (18) 

is the conditional transition probability given that the interarrival time 
6 n = x (constant). For, if f„ = j and d„ = x, then £, 1+] has a Bernoulli 
distribution, either with parameters j + 1 and e~* x when the nth call 
realizes a connection, or with parameters j and e~" z when the nth call 
does not. The system is said to be in state E k at the rath step if f„ = k. 
Starting from the initial distribution {P k (1) } the distributions (P* ) 
can be determined successively by the following formulas: 

Pu n+y) = Z p Jk P/ n) (n- 1,2, ...). (19) 

However, it turns out that in many cases it is more convenient to 
determine the binomial moments of {P k (n) } first. Bj' definition, 



*--.{fr)}-£®*« 



(r = 0,1,2, •••) (20) 
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is the rth binomial moment of {P* (n) }. If we suppose that U, < Ci/r\ 
where d is a constant, then it can be proved that every U T (n) exists 
and Ur n) < C/r\ where C is a constant. Thus the distribution [P k in) \ 
is uniquely determined by {tV'l- We obtain from (20) that 

— * f r \ TT (») 



r=k 



Pk M = E(-ir*(i.)^ (B) (fc- 0,1,2, •••)• (21) 



This is the inversion formula of Jordan. 
It is convenient to use the related quantities 

Fp («) = E |^ pjn | = g ^ p ^u) (r = o, 1, 2, • - • ), (22) 
whence by inversion 

p fc P fc (n) = Z(-D r - fc flV^ > - (23) 



Now we shall prove 

Theorem 1. We have U (n) = 1 (n = 1,2, ■■■) and 

U r u+1) =<p r (Ur M + Vr-i ln) ) (n = 1, 2, • ■ ■ ; r = 1, 2, • • ■ ), (24) 

where <p T = <p(rn). Further 

V r {n) = g (j) (A y -> r )C// n) (r = 0, 1, 2, • • • ), (25) 

where 

± j - r vr = Z(-iy^~ r )pj->- (26) 

Proof. First of all we note that the rth binomial moment of the Ber- 
noulli distribution [Q k ] with parameters n and p, that is, that of 



Q* = (fc)p (i ~p) (fc-0,1, 



5 r = £ (^ <?* = (?") P r (r = 0, 1, • • •, »). (27) 



is given by 

— V7 V. 

Using (27), we get by (18) that 
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whence 

E<l I t - S, 



= (p r 



[-('!■)♦•©] 



(28) 



If we multiply both sides of (28) by P, (n) and add them for every j, 
then we get (24). We obtain (25) if we put (21) into (22). This com- 
pletes the proof of the theorem. 

Starting from U r ' (r = 1, 2, • • •) the binomial moments U, (n = 
2, 3, • • • ) can be obtained recursively by (24) and (25). If, specifically, 
£(0) = i and ri = x then & has a Bernoulli distribution with parameters 
i and e - ' 1 * and thus, for £(0) = i, 



in 



ur = e 



Pr 



(r = 0,1,2, •••)• 



\ r /J V'y 

Remark J. If we introduce the generating functions 



U r (w) = 2 tV 



(n) 



W 



and 



7r(») - Z Vr M vf 



(29) 



(30) 



(31) 



and suppose that £(0) = «, then by (24) and (29) we get that 



U r (w) = 
and evidently 



lVp r 



1 — Wipr 



+ V r -l(w] 



Uo(w) = 



1 — w 



(r = 1,2, •••), (32) 



(33) 



Note also that (21) implies that 



r-t /?' 



W»). 



(34) 



Example 1. In the infinite line case, i.e., when p } = \{j = 0, 1, 2, • ■ • ), 
v u) = jjU) and y T ^ = jj^ for r = o, 1, 2, • • •. If we suppose 

that f (0) = «', then by (32) we get 



1 - u*> r _\'7 J 



(r=l,2,-..) (35) 
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and Uo(w) = wj{\ — w). The .solution of these equations is given by 



U r (w) = n 



W(fij 



\tt/1 — wip> 



En 



y-0,1,2,.-.), 



(>-0 \1 — Wipj/ J {j=0\J/ y~0\ W<P» 

where the empty product means 1. The distribution [P k ln } is deter- 
mined by (34). 

Example 2. For a loss system with m lines, i.e., when p,- = 1 (j < m) 
and pj = (j ^ m), in the case f(0) = i £ m we have 



(«) 



r (n) 



t/ m; __ n <n; 



r„ 



(/• = 0,1, 2, ..,m- 1) 



and 



y w = yw = ( r = m,»i+ 1, •■•)• 



Thus, 



7,(W) = tfr(«) - W ^-(^ ( '' = °> *' 2 ' • ' ■' m " 1} 



and 



F r (w) = C/r(w) = (/• = m, m 4- 1, • • •)• 
By (32) we get 



tfr(w) = 



W<p r 



1 — Wty> r 



t 4- c/ r _i(«0 - ^ r '" 1 ) U m (w) 



(36) 



(r = 1,2, ■••,»») 



and Uu(iv) = w/{\ — w). The solution of these equations for r = 0, 1, 
2, • • • , in is given by 



UAiv) = 



T r (w) 



E , 



LSu/ri(w)JLS\J/r w («)J 



[37) 



where 



r r («) = n 



w;^, 



(r = 0, 1,2, 



,_0 \1 — Wipi/ ' 

and r_i(«») - 1. Finally, (P* 00 ) can be obtained by (34). 
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V. THE LIMITING DISTRIBUTION {P k } 

Iii the Markov chain {£„} the states E , E\ , ■ ■ • , E m form an irre- 
ducible closed set, while E m , E m+ i , • • • are transient states. If either 
m = oo or to < <*>, but we restrict ourselves to the states E , Ei , • • • , 
E,„ , then the Markov chain {£„} is irreducible. The Markov chain {£„} 
is always aperiodic. Accordingly 

lim/V' = P, (ft = 0, 1,2, ..-) 

n-»oo 

always exists and is independent of the initial distribution. There are 
two possibilities: either every P k = (k = 0, 1, 2, ■••) or {P*J is a 
probability distribution. (In the second case P k > if k ^ m and P k = 
if A; > m.) In the second case [Pk] is the unique stationary distribu- 
tion of the Markov chain (£„} and conversely if there exists a stationary 
distribution then it is unique and agrees with the limiting distribution 

{Pk}. 

In the particular case p,- = 1 (j = 0, 1, 2, • ■ • ) the limiting distribu- 
tion always exists, as has been proved in Ref. 10. In this special case 

A= Z(-i) r cv> 0. 

If we consider an arbitrary sequence [p 3 ] then evidently 

^o ^ Z (-i)U > 0, 

r=0 

whence it follows that {£„} belongs to the second class; that is, \P k ] is a 
probability distribution. 

The stationary distribution [P k } is uniquely determined by the fol- 
lowing system of linear equations: 

Pt = E VikPi (38) 

and 

£P*=1. (39) 

k=\ 

Since in this case P k " = P k for every n, we get (38) by (19). Now let 
us introduce the binomial moments 

U r =t(f)Pk (r- 0,1,2,-..) (40) 
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and define 



Vr = £(*)P*P*. 

k=r \T / 

By inversion we get, from (40), 

>r-A- (r 



(41) 



Pk = E(-D (I ^ (* = 0,1,2,.-.) (42) 



and similarly, from (41), 

P* 



p.-Ec-ir'ilK. 



(43) 



The binomial moments t/ T (r = 0, 1, 2, • • • ) can be obtained by the 
following 

Theorem 2. We have U = 1 and 



Ur-z^—Vr-i (r= 1,2,..- 

1 — Vr 



where <pr = <p(rn). Further, 



V r = E (!) (A^Pr) f/y (r = 0, 1, 2, • ■ • ), 



where 



A'- r p, = E(-D' 



./ - r 



P/-» 



(44) 



(45) 



:46) 



Proof. This theorem immediately follows from Theorem 1 if we put 
U r in) = Vr, TV"' = r r in (24) and (25). 

Remark 2. In many cases there is a simple relation between the gen- 
erating functions 



u(z) = E p >J 



(47) 



and 



V(z) = E P*PJ (48) 

A-=0 

when U r (r = 0, 1, • • •) can easily be obtained by (44). For, 

p ipnTUfl ( r = 0,1,2,.-.) (49) 

r ! L dz T _ z=i 
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and 

Theorem 3. The binomial moments U r (r = 0, 1, 2, • ■ •) satisfy the 
following system of linear equations: 

E (-ir* (V) (W f - ^^±- J t/ r+1 ) = (r = 0, 1,2, • • •) (51) 

and 

Ur+1 = ^±1_ £ (A ( A ^ Pr )Uj (r - 0, 1, 2, • • •), (52) 

?y/ie/'e k'~ r p r ts defined by (40). 

Proof. If we put (42) into (43) and use the relation (44) then we 
get (51). If we eliminate V, from (44) and (45) then we get (52). 

Remark 3. If p m = then U r = for r > m, and in this case, start- 
ing from U m , the unknowns U m -i , U m -2 , • ■ • , U can be obtained 
successively either by (51) or by (52) and finally £/ = 1 determines 
U m . If the higher differences of p r vanish, then (52) can be used suc- 
cessfully for the determination of the binomial moments U r . 

Example 3. If p, = 1 0' = 0, 1, 2, ■■■) then V r = U r (f - 0, 1, 2, • • ■ ) 
and, by (44), 



£/, = r -^- t' r _i (r = 1,2,-..), 

1 — ¥»r 



whence 



tfr=Il(V-^-) (r-1,2,-0 (53) 

and Ua = 1. The distribution {P k \ is given by (42). 

Example .',. Let p> = I if j < ??* and p,- = if j ^ m. Then 

K, = U r -(™) U m (r = 0,1, •••,w) 

and 

7 r = u r = (r = m -f 1, m + 2, ■ ■ • ). 



STOCHASTIC PROCESSES WITH BALKING 



807 



By (44; 



and the solution of this equation is 

t ( m ) 1 

v _ c U\i)c, 

Ur — <■ r „, , -. _ 



(r = 1,2, .-•,»), 



S17M 



(/• = 0, 1, • ••, w), 



54) 



where C' r is defined by (3). U T = if r > m. Finally, \P k ) is given by 

Example 5. Let p = 1 and p, = p (j = h%> ' ' ' )» 9i = 3 (j " *> 2 > 
• • • ), where p + </ = 1. Then 

V r = pU r (r = 1,2, •••), 

y = p r + <?Po = 1 - q (U, - Ui + U a - • • ■ )• 

Putting V T into (44) we get 



and 



Ur = 



U r = 



V-Pr 



1 — <fr 



<P\ 



U T -i (r= 1,2, 



- q(Ui - U 2 + U 3 - •••)]■ 



l—^i 
The solution of this system of linear equations is 



v r = 



(r = 0, 1,2, •••), 



1-«I (-p) J Cj 



where C r is defined by (3). Finally, \P k ] is given by (42). 

Example 6. Let po = 1, p, = p, and q, = g if J = 1, 2, 
where p + 7 = 1, and p> = if / > m. Then 



(55) 



m — 1, 



Vo - p + 3 P - p/'„, 

- p + gltfc - r, + r 2 - • • • + (-1 )'"L T m ] - pC/ m , 



in 



V r = pUr - pi J U m (r = 1,2, ■••.m), 
V r = Vr = (r - m + 1, m 4- 2, •••)• 
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Now Uo = 1 and, by (44), 
r »i / \ i "• »• / \ i 

E (") i - «23 (-Dr, P 'X ('") jL < 56 > 

;-o \ J / t >P' /-a ,--> \ I / C ip l 

(r = 1,2, ...,m). 

The distribution {P*} is given by (42). 

Example 7. If, in particular, F(.r) = 1 - e~** for .r ^ 0, then <p(s) = 
X/(X + .9) and <p r = X/(X + rp)(r = 0, 1, 2, • • • ). In this case by (24) 
we have 

rnV, = XF r _, (r = 1,2, •••), 
whence 

M c/'(2) = xr( 2 ). 

Forming the coefficient of z ~ 1 we obtain that 

nkP k = Xp^/V, (k = 1, 2, • • • ), (57) 

whence 






Pfc = /J ° ^r PoPl " " " Pk (k = °' l ' 2 ' ' ' ' } ' 

and Po is determined by the requirement that 

x;p* = i. 

VI. THE TRANSIENT BEHAVIOR OF {£(/)} 

In this section we suppose that £(0) = i always. Denote by M ,{t) 
the expectation of the number of calls occurring in the time interval 
(OVJ which find exactly,/' lines busy. Let 

njis) = f e-UMiit), (58) 

Jo 

which is convergent if (R(s) > 0. Now we shall prove the following 
Lemma 1. Define 

*r(«) = g (fy «,(•) (r - 0, 1, 2, • - •) (59) 



STOCHASTIC PROCESSES WITH BALKING 



N0«) 



*r(«) - Z(«JPiw(«) (r- 0,1,2,-..), 



and 



3=r W 

w/w'c/i are convergent if (R(s) > 0. T/iew 

<*>(«) 



*o(s) = 
and if £(0) = i then 



*r(«J = 



1 - #»(«) 



^>(s + r/x) 



1 - p(s + r/i) 
Proof. Since evidently 

00 

n = l 

we have 



*,« - S (-) *(•> - 1 » {•""■ (*) 



T^r \r 



and .similarly 



i-r V / »— 1 V \ ' / 

Now we shall prove that 



E{e 



,-.r B+I In 



£» = J, On = .C, T n = IJ 



-M't 1 )**©] 



g-»-#« e -»(*+»> 



(60) 



(61 



[(j) +***(«)]. (62) 



(63) 



(64) 



(65) 



This follows from the fact that under the given condition £„+i has a 
Bernoulli distribution either with parameters.;' + 1 and c~" x when the 
nth call gives rise to a connection, or with parameters j and f 1 * when 
the nth call does not. Unconditionally we get 



E{e 



M. 

-„ (. + r„) [E {c "■(«;)} + E {«" »- ( f * J ft.}]. 



(66) 
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If {(0) = i then 



E{e~ aTl 



<p(s + rp). 



If we add (66) for n • = 1,2, • • • and ((57) then we get 



* r (s) = tp(s 4- r/i) 



+ * r (s) + * r _ x (s) 



(67) 



(68) 



(r = 0,1,2, ...), 



where ¥_i(s) = 0. Thus we get (61) and (62). In many cases use of 
Lemma 1 determines $ r (s) (r = 0, 1, 2, • • •) explicitly. 
Remark 4- From (59) we obtain by inversion 



M*(«) = Z(-D 



r-* / r 



*r(s). 



(69) 



The functions n k (s)(k = 0, 1, 2, • • •) can be determined also by the 
following system of linear equations: 



£1> W ' 



> + >>) 



1 - p(S + 77*) _V/ fc-r-1 V " 1. 



I: 



(70) 



which we get if we put (59) and (60) into (62). 

If we know #,(s) (r = 0, 1, 2, • • • ) then P*(0 can be determined by 
the following 

Theorem 4- The La-place transform II* (s) is given by 



n,( s ) = Z(-i) 



r-k / r 



/3 r (s), 



where 



PM = 



[1 - <p(s + r/i)]*r(s) 



(r = 0,1,2, 



<p(s + r/t)(s + r/x) 
Proof. Let the rth binomial moment of JPa-(/)J be denned by 



(71) 



(72) 



B,(0 = E 



a-sC)™ 



(r = 0,1,2, •••)• (73) 



By using the results of Ref. 10 we can see that B T {t) ^ C/rl for every 
t ^ 0, where C is a constant. Thus the probability distribution {Pk{t)\ 
is uniquely determined by its binomial moments. From (73) we get by 
inversion 



,r-k t r 



Pt(t) = E(-D M (£)/>' 



(74) 



STOCHASTIC PROCESSES WITH BALKING 811 

If 

r (s) = I e - at B r (t) (It 

Jo 

and we form the Laplace transform of (74), we get (71). Now let us 
determine r (s) (r = 0, 1, 2, • • • ). 
If $(0) - t, then 

Br (t) = (l\ e- r "[l - F{t)] 

(75) 

+ g [pi ('' + ') + <n (§\ { •" n,M [l " nt - «)] eHfyOO, 

where itf/(t) is denned by (63). For, if there is no call in the time inter- 
val (0,t] then £(/) has a Bernoulli distribution with parameters i and 
e' 11 '. If the last call in the time interval (0,/] occurs at the instant u and 
in that instant the number of busy lines is j, then (■(<) has a Bernoulli 
distribution, either with parameters j + 1 and e - *'""' when this call 
gives rise to a connection or with parameters j and c~' i( ' _ ' ,) when this 
call does not. If we also take into consideration that the last call oc- 
curring in the time interval (0,/] may be the 1st, 2nd, • ••, nth, • • • one, 
then we get (75). Forming the Laplace transform of (74) we get 

„,(.) _!-*•+ W> id 



i + * IV, 



+5[*('t l ) + 'C0. 



M>(«) 



where j*i(«) i y defined by (58). By using the notations (59) and (GO) 
we can write also that 

M- 1 -?;V , {® + * M+ *'- t( ' ) }- (77) 

Taking into consideration the relation (68) we obtain finally 

(s 4- rn) <p(s + rn) 

which was to be proved. 
Example 8. Define 



cm = n ( . v(a t i\ ) (T-o.1,%-) 

,-0 \1 - <p{8 + IH)/ 



(79) 



812 
and 
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C-xOO - i. 



If p . = 1 (j = 0, 1, 2, • • ■ ) and £(0) = i, then * r (s) = <i» r (s) (r = 
0, 1,2, •••) and, by (62), we get 



*M = 



<p(s + rn) 



+ *- (s) ] 



[1 - p(s + rn)] 
where <&-i(s) — 0. The solution of this recurrence formula is 



(r = 0, 1, ..-), (80) 



*r(«) = C r («) § ( •) 777^ ' 



(81) 



where C r (s) is denned by (79). 

Example 9. If p> = 1 whenj < m and p> = when j ^ OT and £(0) = 
i ^ m, then 







¥ r (a) = $ r (s) - ( ) f» m (s) (r = 0, 1, • • -, m) 



and 



# r (s) - $ r (s) - (r = m + 1, m + 2, • • • ). 
(• + 



By (62) 

[1 - ip{s + rti)] _\r > 
for r = 1, 2, • • •, w. The solution of this equation is 

C r 



"(J + «^W "(."i) •*(•>] (82) 



*r(«) = 






™\ 1 (L£v/Ci(a)jLSv/C*-i(8)J 



(83) 



-i(«) 



" LS \j)cXs)\ IS* \j) c~i 

where C/(s) is defined by (79). 

VII. THE LIMITING DISTRIBUTION {Pa-*} 

Now we shall prove 

Theorem 5. If F(x) is not a lattice distribution and its mean a is finite, 
then the limiting distribution 



lim P k (t) = P k ' 

t-»o0 



7c = 0,1, •••) 
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exists and is independent of the initial distribution. We have 

p *+* = Tirrk (*- 0,1,2,-..) (84) 

(k + l)ap 
and 

Po* = 1 - - £ ^ , (85) 

where {P k } is defined by (38). 

Proof. By the theory of Markov chains we can conclude that 

lim*^-* (86) 

<_oo t a 

Furthermore, it is clear that the difference of the number of transitions 
Ek — > Ek+i and E k +i — » 2£* occurring in the time interval (0,t\ is at most 
1. Accordingly, if we denote by N k (t) the expectation of the number of 
transitions E k +i — * E k occurring in the time interval (0,/], then 

\p k M k (t) - N k (t)\ ^ 1 (87) 

for all t ^ 0. Further, 

N k (t) = (fc + l)ju f 2Vi(m) du, (88) 

for, if we consider the process {£(/)! on b' a t those instants when there 
is state E k+ i , then the transitions E k+i — > E k form a Poisson process of 
density (k + 1 )/i. Thus, by (86), (87), and (88), 

.. (k + l)v[* , ,, r N k (t) .. Pk M k (t) p k P k 
lim ■ — — / P*+i(w) dw = Inn — - — = lim- — = - — ; 

(^oo t Jo <-» t t-,*> t a 

that is, 

lim I f P k+i (u) du = ., V _f" (k = 0, 1, 2, • ■ ■ ). (89) 

t-oo t Jo (k + l)an 

If we prove that the limiting distribution 

limP*(0 = P k * (k = 0, 1,2, •••) 

t-»w 

exists, then it follows by (89) that 

*V* = // ZTS (fc = 0,1.2, •••), (90) 

(A; + l)a/i 
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and so 

A *= l - t Pw* - 1 - - tjr^u • (91) 

t=o an k=o {K + I) 

To prove the existence of the limiting distribution we need the follow- 
ing auxiliary theorem: If F(x) is not a lattice distribution, then 



lim 



PkM k (t + h) - PjMkit) (92) 



exists for every h > and is independent of h and the initial state. This 
is a consequence of a theorem of Blackwell. 14 For the time differences be- 
tween successive transitions E k — > E k+i are identically distributed, in- 
dependent, positive random variables, and, if F(x) is not a lattice dis- 
tribution, then these random variables have no lattice distribution 
either. If (92) exists, then it follows that 

{ ^MUt + h) - MM = ]im MUt) = Pj (t _ M , V ..). (93) 

(-00 h t-*oo t a 



Now, by the theorem of total probability, we can write 

(94) 



P k (t) = (*) e- k " x (1 - e'" x y~ k [1 - F(t)] 



+ E f *jk(t - «)[1 - FH ~ u)] dMjiu), 



where v Jk (t) is denned by (18) and it is supposed that £(0) = i. The 
event £(0 = k may occur in several mutually exclusive ways: there is 
no call in the time interval (0,i] and, with the exception of k, all the i 
connections terminate by t; or the last call in the time interval (0,/] is 
the nth (n = 1, 2, • • • ) one and it finds state Ej ij = k — 1, k, • • • ). 
If t„ = u (0 < u ^ t), then during the time interval iu,t] no new call 
arrives [the probability of which is 1 — Fit — u)] and with the excep- 
tion of A; connections every connection terminates by t [the probability 
of which is TTjk(t — u)]. 

Applying Blackwell's theorem to (94) and using a < <x>, it follows 
that 



Km P»(l) = P k * (fc = 0,l, •••) 



t- 



exists and 



i=fc-i 
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where 

Pik* = - f »#(*)[1 - F(x)] dx. (96) 

It is easy to see from (95) that |P*-*) is a probability distribution. 

VIII. THE DETERMINATION OF 7* (s) 

Define 

yM = r e -"dG k (x) = g*^, (97) 

where D (s) = 1. We are going to determine D r (s)(r = 1, 2, • • • ). 
Write D r (s) in the following form: 

DM = zO^Dois), (98) 

where A J Z) (s) is thejth difference of D r (s) at r = 0; that is, 

A J 'A,(s) = £ (-I)'"'' M AC*). (99) 



Then D r (s) is uniquely determined by its differences. 

Now we shall prove 

Theorem 6. Starting from D (s) = A°A>(s) = 1, the functions D r (s) 
(r - 0, 1, 2, • • •) and the differences A'i>o(s) (j = 0, 1, 2, ■ ■ • ) can be 
obtained successively by the recurrence formulas 

±, {-!)-' {*\ DM 

W . '100) 



AU( S ) = = gii + jl £ (^ ( A-p,)A' +1 Do( S ) ( 101 ) 
1 - <p(s + jn) ,=o VV 

respectively. Here 

A J -' Pl — £(-D'f J 7 *W- (102) 
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Proof. By the theorem of total probability we can write for r = 0, 1, 
2, • ■ • that 

G r (x) = ft(.)«"* v (l -«"") H 

•feU^-i/)*-*^-?/) (103) 

+ qfijix - y) *■■■ *G r (x - y)} dF(y), 

where the empty convolution product is equal to 1. Let us consider the 
instant of a transition E T -\ — » E T and measure time from this instant. 
Then G r (x) is the probability that the next transition E r — » E r+i occurs 
in the time interval (0,x\. This event may occur in the following mutually 
exclusive ways: the first call in the time interval (0,x] arrives at the 
instant y(0 < y ^ x), it finds state Ej (j = 0, 1, • • ■, r), the prob- 
ability of which is 



". «"*'(! - c-""Y~ } , 



and, in the time interval (y,x], a transition E r — > E r+ \ occurs, the prob- 
ability of which is 

PjGj+i(x - y) * ■■■ * G r (x - y) + qfifa - y) * • • • * G,(x - y). 

Introduce the notation 

q r j(s) = (\\ f e-'^-^il - e-^Y' 3 ' dF(x) (104) 

and form the Laplacc-Stieltjes transform of (103); then 
Yr(s) - £**<•)[» II y<(s) +q,UyM~\ (r =0,1,2,-.-), 

where the empty product is 1. Now using (97) we find 

DM = £ q r j(s)[ Pj D ]+l (s) + ff/Dy(«)] (r = 0, 1, 2, • • • )• (105) 

This is already a recurrence formula for the determination of D r (s) 
(r = 0, 1, 2, • • • ), but the coefficients can be simplified further. 
If we form 

A'A,(.s) = J2(-l)*" l ((jDM, 
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where Di(s) is replaced by (105), and take into consideration that 

s ( " i)j_ 'C) fy ' ,,(s) = ( - i) "(j)» (a+ * ,) ' (io6) 

then we obtain 

A%(s) = <p(s + jy.) E (-l)^(y W>*w(») + g«D«(«)]. (107) 

Now, comparing (99) and (107), we obtain (100). 
On the other hand, by (107) it follows that 

A J 'Do(s) = <p(s + Jv)A%(s) + *>(« + » E (-D^RJ^W, 

1=0 \i/ 



whence 



and here 



A'DbW = . y(8 +^ l) . . AlpoADo(s)] 

1 - p(s + Jn) 



(108) 



where 



A J '[poADo(.s)] = g Kj (A^>,)A ,+, Do( S ), 



A-v^E(-i)f ;^. 



(109) 



(110) 



This proves (101). 

Example 10. In the infinite line case, i.e., when pj = 1 (j = 0, 1, 2, 
• • • ), (101) has the following simple form: 

A^DoOO = IZJ^tp- A'Z>.(«) (> -0,1,2,.-.), (HI) 



whence 



and 






(112) 



^,)-t(^fl( L^i±^ ). (113) 

;-o V/ j-o \ <p(s + in) J 
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Example 11. If p = 1 and p t = p (j = 1,2, • • • ), then (101) reduces 
to the following difference equation : 

A »< DM - 1 -"<f+f ) A'DM 



+ ( _ l)j g [l-^)] =Q y. 0,1.2,-..). 

A simple calculation shows that the solution of (114) is 



(114) 



A J D (s) = p 



44 L *>(*- 



+ Ui) _ 

q[l-<p(s)] ] ^ ( y fr 1 f l - y (« + fr) l\ 

p*>(«) £i * ; if+i L w(« + to JJ ' 

and finally, 



(115) 



DM = £( . AU( S ). (116) 

Theorem 7. Suppose that £(0) = and under </ws condition denote by 
M k {t) the expectation of the number of calls arriving in the time interval 
(0,t] which find exactly k lines busy. Let 

fi k (s) - [ e~ st dMM. (117) 

Jo 

Then 

pM = 1- n * r-r , (118) 

p fc D fc+1 (s)/x*(s) 

where D, : +i(s) is given by Theorem 6 and nu{s) is given by 

r 

,k 



nM = E(-irM[)*r(s), (119) 



where $ r (s) can be obtained by Lemma 1. 

Proof. The expected number of transitions E k —* Ek+i occurring in 
the time interval (0,/j is evidently p k M k (t). The time differences between 
consecutive transitions E k — > E k+X are identically distributed, independ- 
ent random variables with distribution function R k (x). By using renewal 
theory we can write that 

p k M k (t) = G (l) *G l (t) * ■■■ *G k (t) 

(120) 
*[I{t) + R k (t) + Rfc(f) */2 t (0 + -..], 

where /(/) = 1 if t ^ and I(t) = if / < 0. Forming the Laplace- 
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Stieltjes transform of ( 121 ), we obtain 

1 - Pk {s) D k+l (s)\\ - p k (s)\ 

whence (118) follows. 

Since we know the distribution functions G k (x) and R k (x) (k = 0, 1, 
2, • • • ), the distribution of the number of transitions E k — > E H +\ occur- 
ring in the time interval (0,t] can be obtained easily. 

IX. THE OVERFLOW TRAFFIC 

Suppose that pj I (./ = 0, 1, 2, •••) and that the telephone lines 
are numbered by 1, 2, 3, • ■ • . Further suppose that an incoming call 
realizes a connection through the idle line that has the lowest serial 
number. Consider the group (1, 2, ■ ■ ■, m). Denote by ir m " the proba- 
bility that the nth call finds every line busy in the group (1,2, ■ • •, m). 
The distances between successive calls which find every line busy in the 
group (1, 2, ■ • ■ , m) are evidently identically distributed, independent 
random variables with distribution function, say, G m (x). 

Palm 2 proved that 

,• (n) ■'■ 

lim*. , x (122) 



S(:)i" 



where C T is defined by (3). This is in agreement with (6). In this case 
it is easy to see that 7r„, < " ) = P,,/"', where the distribution {P k ln) \ is de- 
fined in Example 2 of Section IV. 
In Refs. 10 and 11 it is shown that 



.'u 



h \/-/UL <p(s + in) J 

T ( m + ] \fi f 1 - ^ s +v0 l ' 
h \ r yUL *»(* + tM) J 



123) 



where the empty product means 1. It is easy to see that G m (x) agrees 
with the corresponding G m (x) defined in Section VIII when p> = 1 
(J = 0, 1, 2, • ■ ■ ). Thus (123) can be obtained from (07) and (113). 

Remark 5. Denote by T„, the expectation of the random variable 
which is the difference of call numbers of successive calls, both of which 
find all lines busy in the group (1. 2, ■••, m). Knowing T m , we can 
write that 

■w- Km.fi =~ (124) 

B-»°0 1 m 



820 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1961 

and 

f xdG m (x) = aT m . (125) 

Jo 

In Ref. 6 it is shown that r r (r = 1, 2, • • • ) satisfies the following 
recurrence formula: 

r r = ? r . (ri + r 2 + • • • + r r ) + ?r>1 (r 2 + r 3 + • • • + r,) 

(126) 
+ ••• + q r ,r-2(Tr-i + r r ) + ry r . r _,r r + 1, 

where 

q,j = Q j~ r*-(i - rr'm*) u = o, 1, • - ( r). (127) 

The solution of (126) is given by 

r, = ± (') IT i 1 -^) (r-1,8,.-). (128) 
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